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Abstract. Two different "wave chaotic" systems, involving complex eigen- 
values or resonances, can be analyzed using common semiclassical methods. 
In particular, one obtains fractal Weyl upper bounds for the density of reso- 
nances/eigenvalues near the real axis, and a classical dynamical criterion for a 
spectral gap. 

1. Introduction 

These notes present a sketch of semiclassical methods which can be used to 
describe the spectral properties (and as a consequence, the long time properties) of 
a certain class of 1-particle quantum chaotic systems. Here are two examples: 

• damped waves i/j(x, t) on a compact riemannian manifold X of negative sec- 
tional curvature. The dynamics is described by the damped wave equation 

(1) {d 2 t - A x +2b(x)d t )iP(x,t) = 0, 

and the damping function b(x) > is assumed to be smooth. 

• quantum scattering on R d , described by the Schrddinger equation 

(2) imi>(x,t) = p(ttyi>{x,t) , p(h) = -^ + v[x), 

where the potential V(x) has compact support and consists of 3 peaks 
centered on an equilateral triangle; h is Planck's constant. 

These two systems seem very different. The configuration spaces are respectively 
compact and of infinite volume, the wave equation does not depend on K; the 
damped wave equation can be rewritten in terms of a contracting semigroup, while 
the propagator e — , *- p (ft)/fi i s unitary. 

In the damped wave situation, Eq. ([1]) can be diagonalized by a discrete set 
of metastable modes e~ ikjt ipj(x), where tpj(x) G L 2 (X) satisfies the generalized 
eigenvalue equation 

(3) (A + k? + 2ib(x)kj)il;j(x) = 0. 

The eigenvalues kj are complex, and lie in the strip — 2max(6) < 3fc 3 < 0. We 
are interested in the distribution of these eigenvalues in the high-frequency limit 
^tkj — » co. 

In the scattering situation @, the Hamiltonian P(H) is selfadjoint on L 2 (R), with 
absolutely continuous spectrum on R + . However, the Green's function (P(h) — 
z)~ 1 (x,y), well-defined for 3z > 0, admits a meromorphic continuation through 
R + to the lower half-plane, with discrete poles {zj(h)} of finite multiplicities (the 
resonances of P(h)). We will investigate the distribution of these resonances, in 
the vicinity of a fixed energy E > 0, in the semiclassical limit h — > 0. To each 
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resonance is associated a metastable state ipj(x), which is not in L 2 but formally 
decays with time as e~ ltz ^ h ^ h ij)j(x). 

The lifetime of a metastable state is given by the imaginary part of the eigenvalue, 

T = 2p%T> r(3S P- T 3 = 2\<$Sz j {K)\ ■ 

Here are some common features of the two systems. In ij2.ll we show that the 
high-frequency limit 3lkj —¥ oo is similar with the semiclassical limit h — > 0. It 
is then relevant to study the corresponding classical dynamics: in the damped 
wave situation, it is the geodesic flow on X, while in the scattering case it is the 
Hamiltonian flow generated by the Hamiltonian p[x, £) = ^1 — h V{x), in some 
energy interval [E — 8, E + 5]. Our assumptions on X or V(x) imply that these 
classical flows are both "strongly chaotic" . 

In both cases, the long time properties of our quantum system involves a spec- 
trum of complex "eigenvalues" associated with metastable states. Our main aim is 
to understand the distribution of the lifetimes Tj in the semiclassical/high energy 
limit, especially the ones which are not infinitesimally small when h — > 0: we will 
thus focus on resonances such that |£rzj(7i)| = 0(h). 

2. Transformation to nonselfadjoint spectral problems 

Each of these two quantum systems can be recast into a spectral problem for an 
associated nonselfadjoint differential operator on L 2 , with discrete spectrum near 
the real axis. 

2.1. Damped quantum mechanics. Following pQ, let us start from the damped 
wave system. The generalized eigenvalue equation ([3]) for %lkj 3> 1 can be rewrit- 
ten using an effective "Planck's constant" fi « (Sftfcj) -1 , and replacing kj by the 

"energy" Zj = ^f- = 1/2 + 0(h): 

(4) P dw (h)^j = zjfy + o(h 2 ), p dw (h) = - ihb(x) . 

The principal symbol of the operator Pdw{K), Po(x, = |£| 2 /2, is real and generates 
the geodesic flow on X. The skew-adjointness of Pd w (K) only appears in the subprin- 
cipal symbol —ihb(x): the latter does not influence the classical dynamics, but is 
responsible for the decay of probability along the flow. Indeed, for -0o a wavepacket 
microlocalized on p = (x ,£o) e T*X, its evolution ijj{t) = e^ %tPd ^ Ks >l h ^)Q will be 
another wavepacket microlocalized at p t = (x t ,^ t ) — $'(^0), with total probability 
reduced by a finite factor 

(5) mie^.^VH- 

In the limit h — > 0, one can speak of a damped classical dynamics: each point p t 
evolves according to the geodesic flow, and carries a weight which gets reduced by 
the above factor along the flow. 

The horizontal spectral density of Pdw{K) is given by Weyl's law, which (to lowest 
order) does not depend on the damping yQ: for any c > 0, 

(6) USpec(P dw (fi)) n ([1/2 - ch, 1/2 + eft] + iR)} = h- d+1 (cC x + 0(1)) , 

where Cx > only depends on X. On the other hand, we will see that the dis- 
tribution of the imaginary parts Qzj (H) strongly depends on the interplay between 
the geodesic flow and the damping b(x). 
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2.2. Complex scaled scattering Hamiltonian and open dynamics. The scat- 
tering operator P{K) in ([2} can be transformed to a nonsclfadjoint one through the 
complex scaling method 2 . One deforms the configuration space M. d into a complex 
contour Tg = {x + i8f(x)}, where f(x) = for x in a ball B(0,R) containing the 
support of V(x) (the "interaction region"), while f(x) = x for |a;| > 2R. We take 
an angle 8 = Mhlog(l/h), M > fixed. This leads to a "scaled" operator Pg(h), 
which is no more selfadjoint: in the sector —29 < arg(z) < it admits discrete 
eigenvalues, which correspond to the resonances of P{h). 

Our quest for resonances has turned into the spectral study of Pe{h). This 
operator admits the symbol 

p e (x, = p(x, - i6{£, df(xYO + 0{6 2 )\i\ 2 , p(x, £) = !SL + V(x) . 

The operator Pe{K) presents similarities with (Q}: its principal symbol is 
real and generates the Hamiltonian flow on p^ 1 (E), while its imaginary part is of 
higher order h\og(l/h), and generates a damping outside B(0,R). One difference 
with Pdw lies in the strength of this damping: for an initial wavepacket localized 
on a point po e p^ 1 {E), the full probability after time t will be reduced by a factor 
|j 2M / K>. 1, /(i.) S«)ds. probability is semiclassically strongly suppressed as soon 
as the trajectory enters the zone where (£, df(x) t £ i ) > (this "absorbing zone" 
contains the exterior of B(0, 2R)). 

Classically, this corresponds to an "open dynamics" : the point pt evolves ac- 
cording to the Hamiltonian flow, but it gets "absorbed", or "killed" as soon as it 
enters the absorbing zone. 

Trapped set. For any energy E > 0, the forward (resp. backward) trapped set K% 
(resp. Kg) is defined as the set of initial points which remain bounded for all 
positive (resp. negative) times: 

Kl = {p Q = (x ,Z Q ) ep-\E), \x t \<R, Vf^O}, 

while the trapped set is made of their intersection Ke = Kg n Kg. Notice that 
Ke is a compact flow-invariant set. The above remark shows that any wavepacket 
localized on a point pa ^ Kg will be absorbed after a finite time, namely the time it 
takes to enter the absorbing zone. On the other hand, a point p € Kg will converge 
to the trapped set Ke as t — > oo. This argument shows that, in some sense, long 
time quantum mechanics (at energy w E) takes place on Ke- 

3. Fractal Weyl laws 

The above argument can be made precise when estimating the number of reso- 
nances of P{K) near E. In the case we are interested in, Ke is a hyperbolic repeller 
(that is, there is no fixed point on Ke, and all trajectories are hyperbolic), this 
number is bounded above by a fractal Weyl law directly related with the geometry 
of K E 0]. 

Theorem 1. Assume that the trapped set Ke at some energy E > is a hyperbolic 
repeller, and write its Minkowski dimension diuiM (K e) = 1 + 2v. Then for any 
c, a > ; one has in the semiclassical limit 

(7) |({Res(P(ft)) n ([E -ch,E + ch] - ih[0, a})} = 0{h- u -°) . 
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This theorem is proved by conjugating Pg(h) by a suitable "weight" G(h), so that 
the symbol pg^c of the conjugated operator e~ G ^ Pge G ^ satisfies 9pe,G < ~2Ch 
outside the \FK- neighbour hood of Ke- Estimating the volume of this neighbour- 
hood, and some involved pseudodifferential calculus on ps,Gi lead to the above 
upper bound. 

A similar argument can be used to study the distribution of decay rates Qzj/h 
for the operator (|4]). For any time T > 0, one can construct a weight Gt{K), such 
that the conjugate operator Pd w ,G T (k) — e~ GT ^ Pd w (h)e GT ^ admits the symbol 

(8) p dw ,a T M = ^-ihb T ( x ,o + o(h 2 ) 1 

where the subprincipal symbol &r(p) = T^ 1 f^^ 2 b(xt) dt is the average of the 
damping along the flow. The geometric assumption of negative curvature implies 
that the geodesic flow on X is Anosov, in particular it is ergodic. This implies that, 
on the energy shell p~ 1 (l/2), the time average br{p) converges almost everywhere 
to the microcanonical average b = Vol(X) -1 J x b(x)dx when T — > oo. From there 
one can deduce that most of the eigenvalues G [1/2, ch, 1/2 + ch] concentrate 

near the "typical line" = — hb [I]. 

Using finer properties of the Anosov geodesic flow, one can estimate the number 
of eigenvalues away from this "typical line" j3]. To state the result, we need to intro- 
duce the extremal ergodic averages of the damping, b— — lim^—^oo minp-1^/2) ^Ti 
and similarly for b + . 

Theorem 2. Assume X is a compact surface of negative curvature. Then, there 
exists a function H : R — > R, strictly concave on [6_, 6+] and equal to —00 outside, 
with maximum H(b) = d — I, such that for any c > and any a G [0, b], 

(9) i{^c{P dw {h)) n ([1/2 - ch, 1/2 + ch] - ih[o, a})} = o{h- H ^-°) . 

A similar estimate holds for the range — 3z > ha, a > b. 

Comparing this "fractal Weyl upper bound" with the Weyl law © confirms that 
most resonances are on the typical line. The above theorem is obtained by studying 
the large deviations of the value distribution of bx in the limit T — > 00: roughly 
speaking, for a < b, the volume of the points p £ p~ 1 (l/2) such that &t(p) < ot 
decays like e 7 ^^' - ^ -1 ^. This volume estimate is then used to get ([9]). 

In both situations, the upper bounds on counting resonances/eigenvalues were 
obtained by deforming the operator Pg (resp. Pdw) by an appropriate microlocal 
weight, and studying the imaginary part of the resulting operator by phase space 
volume arguments, where the classical dynamics plays a prominent role. 

3.1. Spectral gaps. We now present a complementary type of spectral informa- 
tion, which can be obtained by a similar method for these two systems. Namely, 
we want to understand if the lifetimes Tj(H) can be arbitrarily large in the semi- 
classical limit; or on the opposite, if there exists a gap of size oc h between the real 
axis and the eigenvalues/resonances. The presence of such a gap has important 
consequences on the long time properties of the system. 

Let us start with the scattering problem. The gap question can be rephrased 
as: "Are the metastable states able to concentrate on Ke when h —> 0?" The 
answer will result from a competition between, one one side, the fast dispersion of 
wavepackets due to the hyperbolic classical flow, on the other side the "thickness" 
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of the trapped set allowing the state to reconstruct itself through constructive 
interferences. 

A dynamical quantity reflecting this competition is of statistical nature, it is a 
topological pressure associated with the flow on Ke- 

V{-<p u /2) = lim i log ]T exp ( - / ' <p u (pt)/2 dt) . 

Each 7 is a periodic orbit in Ke with period T 7 , and f u (p) is the unstable Jacobian 
of the flow. The above mentioned competition lies in the fact that each exponential 
weight gets very small when T 7 — > oo, while the number of terms grows exponen- 
tially. 

The following gap criterion was first obtained [5] for the case of hard obstacles, 
and then generalized in [6] to smooth potentials. 

Theorem 3. Assume the trapped set Ke is a hyperbolic repeller. If the pressure 

V(— tpu/2) < 0, then for any c, e > and any small enough h, the strip 

[E — ch, E + ch] + ih\P{— ip u /2) + e, 0] does not contain any resonance of P{h). 

In dimension d = 2, the condition T'(—(p u /2) < is equivalent with a purely 
geometric statement, namely the fact that the Hausdorff dimension dim# (Ke) < 2 
(notice that dim/fp _1 (i?) = 3). 

In the case of damped waves, the gap question is nontrivial if fe_ = 0, that is, if 
there exists a flow-invariant subset of p _1 (l/2) with no damping. A result similar 
to the one above was obtained in [7j. In this case, the local decay of probability is 
due to both hyperbolic dispersion and damping. 

Theorem 4. Assume that X has negative curvature, and that the topological pres- 
sure V(—(p u /2 — b) < 0. Then, for any c, e > 0, and h > small enough, the strip 
[1/2 — ch, 1/2 + ch] +ih[P(— b— (fu/2) + e, 0] does not contain eigenvalues ofPd w (h). 

4. Open questions 

Most of the above results are upper or lower bounds. The natural question is: 
"Are these bounds sharp?" The fractal Weyl bound J?]) is conjectured to be sharp 
for a > large enough, a fact which has been tested numerically on a number 
of examples, but could be proved only for a very specific toy model [5]. On the 
opposite, the bounds Q for eigenmodes of the damped wave equation are not 
expected to be sharp for all values of a. The size of the gap itself is believed to 
be larger than the topological pressure bound we gave above. Such an "extra gap" 
was proved for the 3-disk scattering, using advanced estimates on classical mixing 

m- 
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